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ABSTRACT 

We present multi-epoch spectroscopy of “heartbeat stars,” eccentric binaries with dynamic tidal dis¬ 
tortions and tidally induced pulsations originally discovered with the Kepler satellite. Optical spectra 
of six known heartbeat stars using the Wyoming Infrared Observatory 2.3 m telescope allow measure¬ 
ment of stellar effective temperatures and radial velocities from which we determine orbital parameters 
including the periods, eccentricities, approximate mass ratios, and component masses. These spec¬ 
troscopic solutions confirm that the stars are members of eccentric binary systems with eccentricities 
e > 0.34 and periods P = 7-20 days, strengthening conclusions from prior works which utilized 
purely photometric methods. Heartbeat stars in this sample have A- or F-type primary components. 
Constraints on orbital inclinations indicate that four of the six systems have minimum mass ratios 
q = 0.3-0.5, implying that most secondaries are probable M dwarfs or earlier. One system is an eclips¬ 
ing, double-lined spectroscopic binary with roughly equal-mass mid-A components {q = 0.95), while 
another shows double-lined behavior only near periastron, indicating that the FOV primary has a GIV 
secondary (g = 0.65). This work constitutes the first measurements of the masses of secondaries in a 
statistical sample of heartbeat stars. The good agreement between our spectroscopic orbital elements 
and those derived using a photometric model support the idea that photometric data are sufficient to 
derive reliable orbital parameters for heartbeat stars. 

Keywords: Techniques: radial velocities — (Stars:) binaries: general — (Stars:) binaries: spectro¬ 
scopic — (Stars:) binaries: close — Stars: kinematics and dynamics 


1. INTRODUCTION 

In addition to being a planet hunter, the Kepler satel¬ 
lite has uncovered a host of unexpected stellar phenom¬ 
ena, including new and exotic classes of binary stars. 
Most stars form in binary stellar systems (Duquennoy & 
Mayor 1991), making it imperative to understand their 
statistical properties since they expose the physical pro¬ 
cesses inherent in star formation and evolution. Kepler 
data have revealed a new class of eccentric binaries in 
which the primary star undergoes tidal distortions and 
pulsations in response to perturbations from the sec¬ 
ondary star. An early example, KOI-54, was reported 
by Welsh et al. (2011) who observed a light curve that 
brightened by 7% every 41.8 days and exhibited addi¬ 
tional modulations between maxima. They attributed 
the primary brightening to the tidal interaction of a bi¬ 
nary system near periastron, while the higher frequency 
variations were interpreted as tidally induced oscilla¬ 
tions. These objects had been predicted by Witte & 
Savonije (1999), Witte & Savonije (2002), and Willems 
& Aerts (2002), but Kepler allowed for an in-depth ob¬ 
servational analysis of the tidal distortions. Later works, 
such as Fuller & Lai (2012), Burkart et al. (2012), and 
O’Leary & Burkart (2014), have shown that the domi¬ 
nant brightening near periastron passage is likely due to 
resonantly locked m = 2 modes that trace the period. 

Thompson et al. (2012, hereafter TH12) reported the 
discovery of 17 additional such binary systems from the 
Kepler archive; the tally is now over 130 (Hambleton et 
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al. 2013). The multi-modal nature of the light curves 
resembles an electrocardiogram, earning them the name 
“heartbeat stars”. TH12 determined an orbital period 
from the light curve and used the model of tidal distor¬ 
tions presented in Kumar et al. (1995) to ht the eccen¬ 
tricity, angle of periastron, and orbital inclination to the 
shape of each Kepler light curve. While the models of 
TH12 appear to provide good hts to the data, their lim¬ 
ited number of radial velocity measurements precluded 
a definitive measurement of the orbital parameters and 
could not address the velocity amplitudes and masses of 
the stellar components. 

In this paper, we present a three-year campaign of op¬ 
tical spectroscopy for a small sample of heartbeat stars. 
We chose a sample of six stars from TH12 as a pilot 
survey to directly measure orbital solutions and test the 
accuracy of the photometry-based model hts to the or¬ 
bital parameters. We conhrm the binary nature of these 
systems and hnd them to be short-period, eccentric bi¬ 
naries. Thus, we provide an independent confirmation 
of the hypothesis that the unusual light curves of these 
systems stem from tidal interactions of eccentric binary 
star systems. Our six targets were chosen to be repre¬ 
sentative of the sample of 17 systems presented in TH12, 
spanning the range of temperatures and periods while 
being bright enough to obtain high-signal-to-noise obser¬ 
vations with the Wyoming Infrared Observatory 2.3 m 
telescope. Section 2 discusses our spectroscopic observa¬ 
tions and data reduction methods. Section 3 contains an 
analysis of radial velocities from the extracted spectra 
and a discussion of our independent temperature mea¬ 
surements. Section 4 presents best-ht orbital solutions 
for the six systems and a comparison to the solutions 
from TH12. All velocities reported in this work are in 
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Table 1 

Targets 


Name 

R. A. (J2000) 

Dec. (J2000) 

Kp (mag) 

# Obs. 

KIC 3230227 

19 20 27.025 

+38 23 59.46 

9.002 

23 

KIC 3749404 

19 28 19.089 

+38 50 13.60 

10.575 

18 

KIC 4248941 

19 10 24.766 

+39 23 51.63 

12.159 

21 

KIC 8719324 

20 04 27.423 

+44 49 32.87 

11.614 

24 

KIC 9899216 

19 40 38.801 

+46 45 02.23 

10.886 

16 

KIC 11494130 

18 56 06.426 

+49 24 55.50 

10.988 

19 


the Heliocentric frame of reference. 


2. SPECTROSCOPIC OBSERVATIONS AND REDUCTIONS 

Spectra were obtained at the Wyoming Infrared Ob¬ 
servatory (WIRO) 2.3 m telescope over the period 2012 
March through 2014 September and include 49 nights of 
observations with an average of 20 observations per star. 
Table 1 provides a list of the targets, the right ascen¬ 
sion and declination of each object, the Kp magnitude 
(referring to the broad optical Kepler bandpass, which 
corresponds roughly to the Cousins Re magnitude), and 
the number of nights that each target was observed. The 
stars have Kp magnitudes between 9.0 and 12.2. 

All data were obtained using the WIRO-Longslit Spec¬ 
trograph with an e2v 2048 x 2048 CCD. A 2000 line 
per millimeter grating in first order provided a spectral 
resolution of about 1.2 A FWHM in the center of the 
5400-6700 A spectral range. A 1".2 x 110" slit oriented 
at a position angle of 0° yielded a reciprocal dispersion 
of 0.63 A per pixel and resolution R^dSOO at 5800 A. 
Exposures ranged from 600 s to 1800 s in 200-300 s in¬ 
crements, depending on the target brightness and nightly 
conditions. Seeing ranged from about 1".2 to 3" FWHM. 

Reductions followed standard techniques for longslit 
spectroscopy. Spectra were flat-fielded from dome quartz 
lamp exposures and wavelength calibrated using copper- 
argon arc lamp exposures to a rms of 0.02 A (1.0 km s“^ 
at 5800 A). Exposures for a given star on a given night, 
taken sequentially in the observing sequence, were com¬ 
bined to yield one spectrum with signal-to-noise ratio 
(SNR) of over 70 per 0.63 A pixel. All spectra were 
Doppler corrected to the Heliocentric reference frame. 
The mean uncertainty on the night-to-night velocity zero 
point was determined to be 4.1 km s“^ rms based on ex¬ 
tensive observations of OB stars within the Cygnus OB2 
Association obtained at WIRO over the same time pe¬ 
riod (Kobulnicky et al. 2014).^ This uncertainty was 
added in quadrature to the radial velocity uncertainty 
determined during cross correlation, discussed in Sec¬ 
tion 3.1. For the system with the largest velocity am¬ 
plitude, KIC 3230227, we expect the radial velocity to 
change by roughly 6 km s“^ per hour near periastron. 
As the time for a series of exposures for this system is 
10-15 minutes, we would expect the total change in ra¬ 
dial velocity to be less than 3 km s“^. Because this is 
significantly less than our other velocity uncertainties, 
we consider it negligible. 

Figure 1 shows a continuum-normalized spectrum of 
KIC 8719324 (Teff =7600 K) as a representative example 

® This uncertainty is a conservative estimate of the velocity 
uncertainties introduced by temperature changes and flexure in 
the Cassegrain spectrograph. 
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Figure 1. Spectrum of KIC 8719324 (Teff = 7597 K)as a repre¬ 
sentative example of a heartbeat star. Labels indicate the strongest 
features in the spectrum, including Ho, stellar and interstellar 
Na I D AA5889,5895, Si II A6347 A, and Si II A6371 A. Horizon¬ 
tal bars mark telluric O 2 and H 2 O bands. Many other unlabeled 
metal lines are visible in the spectrum, providing significant cross¬ 
correlation power for determining radial velocities. 

of a heartbeat star. The strongest features in the spec¬ 
trum are Ha and Na I D, which contains both stellar and 
interstellar components. Additional weak features from 
low-ionization metals are abundant in the spectrum and 
not individually labeled, with the exception of Si II at 
A6347 and A6371 A. No telluric corrections have been 
applied, so atmospheric features are present in the spec¬ 
tra, principally the O 2 band around 6300 A and the H 2 O 
band around 5900 A. 

3. DATA ANALYSIS 
3.1. Cross correlation Methods 

We used the IRAF task XCSAO"^ to measure radial 
velocities from the spectra using cross-correlation tech¬ 
niques. The XCSAO task allows the user to specify a 
template spectrum against which to perform the cross¬ 
correlation. Initially, we tried to use the highest signal- 
to-noise spectrum for each star in our sample as the tem¬ 
plate spectrum. However, even our best spectra had a 
signal-to-noise ratio just over 100 so the observational 
velocity uncertainties were on the order of 1-10 km s“^. 
As some of these systems, such as KIC 9899216, have ve¬ 
locity amplitudes of ^10 km s“^, such large uncertainties 
preclude secure orbital solutions. We then tried perform¬ 
ing the cross-correlation against synthetic spectra from 
the MARCS (Gustafsson et al. 2008) model stellar atmo¬ 
spheres of the nearest available temperature and gravity 
to those reported in TH12.^ We corrected the vacuum 
wavelengths of the synthetic spectra to wavelengths in air 
and smoothed them by 1.1 A to match the resolution of 
our spectra. To ensure that resulting velocities were not 
sensitive to a specihc template, we performed a test on 

^ XCSAO is a radial velocity analysis task from the Smithsonian 
Astrophysical Observatory package RVSAO. 

® KIC 3749404 was correlated against a template of 7500 K 

with logg = 3.5; KIC 4248941 against a template of 6750 K 

with logg = 4.5; KIC 8719324 against a template of 7500 K 

with logg = 3.5; KIC 9899216 against a template of 7500 K with 

logg = 3.5; KIC 11494130 against a template of 6750 K with 
logg = 4.5. 
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target stars at 6750 K, 7500 K, and 8000 K in which we 
cross-correlated against templates mismatched in tem¬ 
perature by several hundred to a few thousand Kelvin. 
The resulting velocities changed by < 0.2 km s“^ and 
the uncertainties increased by < 0.4 km s“^. Similarly, 
we rotationally broadened a template by 100 km s“^ and 
found that the average change in the measured velocity 
was less than 1.5 km s“^, much less than other velocity 
uncertainties inherent in our spectra. We conclude that 
the derived radial velocities are not sensitive to small 
mismatches in temperature, surface gravity, or rotational 
broadening between the star and template. With these 
synthetic spectra, we can both estimate a systemic veloc¬ 
ity and obtain average velocity precisions of < 2 km s“^. 
The accuracy of the velocities is ultimately limited by 
the 4.1 km s“^ night-to-night uncertainties on the veloc¬ 
ity zero point. 

When correlating the stellar velocities with XCSAO, 
the user can specify specific regions of the spectrum to 
be ignored in the correlation. It is particularly impor¬ 
tant, for example, to exclude the Na I D region because 
it contains both stellar and interstellar components. We 
considered four different combinations of spectral cover¬ 
age in the cross-correlation. Our first cross-correlation 
method did not exclude any lines in the spectrum and 
was deemed to yield an unreliable velocity because the of 
the presence of the interstellar component of the Na I D 
lines. The second method used only the Ha (6530- 
6585 A) portion of the spectrum and was also discarded 
as sub-optimal because of the broad wings of Ha leading 
to large velocity uncertainties. Our third method re¬ 
moved both the sodium (5883-5903 A) and Ha features 
from the spectrum; the errors on the resulting velocities 
became larger in the absence of the strong Ha line. Our 
adopted method for measuring velocities removed only 
the Na I D region from the spectrum and used the cross¬ 
correlation of the plethora of metal lines and the strong 
Ha. The system KIC 3749404 shows some line splitting 
around periastron passage, so a measurement of the sec¬ 
ond peak in the cross correlation gave the velocity of the 
secondary when splitting was evident. For the SB2 sys¬ 
tem KIC 3230227, whose components are nearly equal 
and therefore difficult to unentangle in the cross correla¬ 
tion spectrum, we instead measured the velocities of both 
components using a double Gaussian fit to the Si II A6347 
line. The width and depth of the Gaussian fits were con¬ 
strained using unblended spectra obtained near quadra¬ 
ture phase. The velocities were not measured using the 
cross-correlation described above because of the difficulty 
separating the two peaks of the cross-correlation function 
when the components were partially blended. Typical ve¬ 
locity uncertainties using this approach were 3-7 km s“^, 
comparable to or larger than the night-to-night velocity 
zero point uncertainties. 

Tables 2 and 3 provide the ephemerides for the four 
SBl systems KIC 4248941, KIC 8719324, KIC 9899216, 
and KIC 11494130, including the Heliocentric Julian date 
(HJD—2,400,000) of observation in column 1, the phase 
of the orbit in column 2, measured radial velocities (I4i) 
with uncertainties for the primary in column 3, and 
the difference between observed and computed velocities 
from the orbital solution {0\ — Ci) for the primary in 
column 4. Table 4 provides the provides the same infor- 


Table 2 

Ephemerides for Single-Lined Binary Systems 


HJD-2,400,000 b 

(km s“l) 

Oi - Cl 
(km s“l) 


KIC 4248941 

56011.668 

0.32 

7.8(4.8) 

-3.7 

56015.633 

0.78 

-27.9(4.5) 

-1.1 

56097.702 

0.27 

18.7(4.2) 

-0.9 

56098.721 

0.39 

4.7(4.2) 

3.6 

56103.801 

0.98 

24.4(4.2) 

-2.9 

56107.761 

0.44 

-4.0(4.2) 

0.9 

56118.679 

0.70 

-28.4(4.2) 

-1.6 

56128.816 

0.87 

-8.3(4.3) 

7.2 

56169.687 

0.60 

-16.6(4.3) 

4.5 

56178.756 

0.65 

-18.9(4.2) 

5.6 

56179.653 

0.75 

-33.0(4.2) 

-5.5 

56435.801 

0.38 

-2.3(4.2) 

-4.4 

56436.688 

0.49 

-12.1(4.3) 

-1.6 

56437.682 

0.60 

-21.3(4.3) 

-0.1 

56438.766 

0.73 

-31.2(4.2) 

-3.7 

56445.718 

0.53 

-14.0(4.3) 

1.1 

56447.720 

0.76 

-30.2(4.2) 

-2.8 

56796.706 

0.13 

46.7(4.3) 

-0.3 

56798.739 

0.37 

4.6(4.3) 

0.1 

56799.812 

0.49 

-8.4(4.2) 

2.6 

56804.848 

0.07 

55.4(4.2) 

2.4 

KIC 8719324 

56011.727 

0.53 

-28.5(12.5) 

-2.2 

56015.707 

0.92 

-15.5(4.5) 

2.7 

56097.736 

0.93 

-15.1(4.3) 

-2.0 

56098.747 

0.03 

37.1(4.3) 

-14.9 

56103.825 

0.52 

-32.3(4.3) 

-6.0 

56107.765 

0.91 

-19.7(4.3) 

0.0 

56118.702 

0.98 

18.0(4.3) 

-7.7 

56128.760 

0.96 

18.0(5.9) 

11.1 

56169.713 

0.96 

3.9(4.5) 

-4.2 

56178.780 

0.85 

-16.2(4.3) 

12.3 

56179.640 

0.93 

-13.2(4.3) 

-1.1 

56438.751 

0.25 

-10.4(4.3) 

-0.8 

56440.728 

0.44 

-23.9(4.3) 

-1.1 

56441.769 

0.54 

-31.2(4.3) 

-4.3 

56444.731 

0.83 

-35.5(4.3) 

-6.1 

56445.731 

0.93 

-9.5(4.3) 

3.3 

56447.748 

0.13 

14.6(4.3) 

5.0 

56464.690 

0.78 

-34.8(4.5) 

-4.0 

56468.933 

0.20 

-4.8(4.8) 

-1.3 

56795.870 

0.14 

10.5(4.5) 

4.0 

56796.737 

0.23 

-5.7(4.8) 

1.4 

56798.805 

0.43 

-15.0(5.0) 

7.2 

56799.801 

0.52 

-27.2(4.3) 

-1.0 

56804.792 

0.01 

69.5(4.3) 

13.7 


mation as Table 2 for the SB2 systems KIC 3749404 and 
KIC 3230227, with the addition of Vr 2 and O 2 — C 2 for 
the secondary star in columns 5 and 6, respectively. 

3.2. Temperatures 

To measure temperatures for the six systems, we used 
the MARCS model atmospheres (Gustafsson et al. 2008) 
and smoothed the continuum-normalized model spectra 
to the resolution of the observed spectra. We measured 
the equivalent widths of the Si II A637I.4 and Ha lines for 
both observed and model spectra. The EW of the silicon 
line was measured by fitting it with a Gaussian profile, 
while the EW of Ha was measured by fitting it with a su¬ 
perposition of two Lorentzian profiles. We measured the 
equivalent widths (EW) for the range of MARGS model 
atmospheres between temperatures of 6000 K and 8000 K 
and logg between 3.0 and 5.0. Figure 2 plots the stel- 
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Table 3 

Ephemerides for Single-Lined Binary Systems 


Table 4 

Ephemerides for Double-Lined Binary Systems 


Vrl Ol - Cl 

HJD—2,400,000 0 (km s“^) (km s“^) 


KIC 9899216 


56439.692 

0.93 

5.2(4.8) 

-1.9 

56440.715 

0.01 

6 .6(4.8) 

-0.2 

56442.727 

0.18 

-1.0(4.5) 

3.9 

56446.698 

0.51 

-1.1(5.3) 

2.2 

56461.671 

0.76 

4.3(4.8) 

4.3 

56462.679 

0.85 

3.3(5.0) 

0.9 

56795.801 

0.64 

0.4(9.1) 

2.3 

56796.688 

0.71 

-2.1(5.3) 

-1.1 

56798.772 

0.88 

1.0(4.8) 

-3.1 

56799.701 

0.96 

11.6(5.0) 

1.1 

56804.774 

0.39 

-9.4(5.0) 

-5.2 

56819.682 

0.63 

-1.0(4.3) 

1.1 

56821.706 

0.80 

0.4(4.8) 

-0.5 

56822.696 

0.88 

5.6(4.8) 

1.8 

56825.688 

0.13 

-5.4(4.3) 

-0.9 

56829.697 

0.46 

-6.3(4.3) 

-2.7 


KIC 11494130 


56097.661 

0.88 

-39.4(4.2) 

-2.6 

56098.677 

0.93 

-40.1(4.2) 

-0.2 

56107.735 

0.41 

1.9(4.2) 

0.3 

56118.653 

0.99 

-29.7(4.2) 

-1.2 

56128.706 

0.52 

-1.8(4.3) 

3.1 

56169.727 

0.68 

-18.6(4.3) 

-2.1 

56178.792 

0.16 

-1.8(4.2) 

15.8 

56179.625 

0.20 

11.9(4.2) 

-0.8 

56436.675 

0.75 

-22.7(4.2) 

-0.1 

56438.738 

0.86 

-33.7(4.2) 

0.5 

56439.678 

0.91 

-39.4(4.2) 

-0.3 

56440.701 

0.96 

-35.1(4.2) 

2.1 

56444.686 

0.17 

12.5(4.2) 

-1.2 

56445.688 

0.22 

11.7(4.2) 

-0.1 

56447.763 

0.33 

3.7(4.2) 

-2.4 

56796.859 

0.73 

-20.0(4.2) 

1.1 

56798.789 

0.83 

-30.5(4.2) 

1.1 

56799.721 

0.88 

-36.6(4.2) 

0.4 

56804.759 

0.15 

16.6(4.2) 

2.5 


lar effective temperature versus the ratio of equivalent 
widths EWho / EWsiii for both the MARCS models 
and our data. The gray tracks illustrate the relation¬ 
ship between temperature and EW ratio for five differ¬ 
ent model surface gravities, as labeled, with the black 
track showing the nominal main sequence in this tem¬ 
perature range.® We also plot a quartic extrapolation 
of the last several points on the main sequence track to 
provide an estimate of EW ratios for hotter stars. Each 
vertical line marks the equivalent width ratio of a star 
in our sample. Symbols with error bars on the corre¬ 
sponding vertical line lie at the intersection of the ver¬ 
tical line with the main sequence. The uncertainties on 
the equivalent width ratio are measured from the error in 
the profile fit to the spectral lines; uncertainties on tem¬ 
perature are derived from the uncertainties of the EW 
ratio and the local slope of the adopted main-sequence 
track. Although our data are not of sufficient resolution 
and wavelength coverage to calculate a log g for each star, 
other measurements, such as the photometric measure- 

® We adopt the relation between spectral type and log g recom¬ 
mended by Castelli & Kurucz (2004). For late-A through late-F 
stars typical of our sample, log g is ~4.0. 




Vrl 

Ol - Cl 


O2- 

-Cz 

HJD-2,400,000 

4> 

(km s“l) 

(km s“l) 

(km s“l) 

(km 

s-i) 


KIC 3749404 

56439.706 

0.32 

15.4(4.3) 

-0.8 




56440.687 

0.37 

17.2(4.3) 

-0.9 




56442.709 

0.47 

-5.3(4.3) 

-0.6 




56446.711 

0.66 

-42.3(4.3) 

-1.3 




56458.683 

0.25 

8 .4(4.3) 

-3.9 




56460.665 

0.35 

12.9(4.3) 

-4.7 




56795.770 

0.84 

-16.3(4.5) 

1.6 




56796.846 

0.89 

-9.4(4.2) 

3.8 




56798.826 

0.99 

-5.1(4.3) 

0.4 




56799.734 

0.03 

-1.9(4.3) 

0.3 




56804.870 

0.29 

17.6(4.3) 

3.3 




56871.663 

0.57 

-69.3(4.3) 

-5.1 

67.8(35.8) 


-0.3 

56876.656 

0.82 

-17.3(4.2) 

2.6 




56903.726 

0.15 

4.4(4.3) 

-1.2 




56907.638 

0.34 

18.0(4.3) 

0.7 




56913.636 

0.64 

-44.9(4.3) 

1.4 




56916.608 

0.78 

-19.7(4.2) 

3.9 




56931.791 

0.53 

-86.1(4.3) 

-3.1 

95.3(8.7) 


1.4 

KIC 3230227 

56011.659 

0.44 

-16.5(7.7) 

-5.7 

-19.0(25.0) 


1.8 

56015.601 

1.00 

-90.5(5.9) 

-12.2 

46.9(7.7) 


-4.2 

56097.685 

0.64 

19.7(4.5) 

-4.1 

-47.8(5.0) 


9.9 

56098.703 

0.79 

44.3(4.8) 

-3.4 

-83.6(5.6) 


-0.5 

56103.792 

0.51 

-7.3(7.1) 

-8.4 

-17.4(10.0) 


16.2 

56107.717 

0.06 

-146.5(5.0) 

-13.7 

122.6(6.3) 


13.4 

56118.671 

0.62 

6 .9(5.3) 

-13.0 

-48.1(6.3) 


5.4 

56128.722 

0.04 

-116.6(6.7) 

21.3 

94.8(8.3) 


-19.7 

56169.673 

0.85 

58.8(5.3) 

2.0 

-92.6(6.3) 


0.3 

56178.744 

0.14 

-97.5(5.3) 

-2.0 

95.7(6.7) 


26.2 

56179.619 

0.26 

-42.0(5.3) 

9.1 

36.8(6.3) 


14.7 

56435.788 

0.59 

25.9(5.0) 

10.6 

-40.5(6.3) 


8.2 

56436.711 

0.72 

30.6(5.6) 

-6.5 

-51.8(7.1) 


20.1 

56437.705 

0.86 

61.7(4.5) 

3.7 

-84.7(5.0) 


9.5 

56438.716 

0.01 

-99.0(6.3) 

-6.0 

61.0(8.3) 


-5.7 

56439.722 

0.15 

-98.7(5.0) 

-8.3 

62.3(5.9) 


-1.7 

56443.716 

0.71 

9.1(5.3) 

-26.9 

-100.0(6.3) 


-29.3 

56444.714 

0.86 

52.3(5.9) 

-5.1 

-95.8(7.1) 


-2.2 

56445.709 

1.00 

-70.1(4.8) 

-0.4 

43.6(5.6) 


1.7 

56446.952 

0.17 

-81.4(5.3) 

-1.6 

46.6(6.3) 


-6.1 

56447.698 

0.28 

-46.0(5.3) 

0.0 

7.8(6.3) 


-9.0 

56467.678 

0.11 

-110.1(5.0) 

-2.8 

83.0(5.9) 


0.9 

56470.708 

0.54 

34.4(5.6) 

27.1 

-32.6(4.8) 


7.4 


ments in the Kepler Mikulski Archive for Space Tele¬ 
scopes (MAST) archive, suggest that these stars have 
log 5 = 4.0 ± 0 . 2 . As our grid of log 5 values is much 
coarser than this, we assume a main sequence luminos¬ 
ity class for these objects. While this assumption may 
have implications on the temperature, it will have min¬ 
imal impact on the mass. Thus, we advocate that our 
temperatures be used with caution, but find them suffi¬ 
cient for the present analysis. Most stars fall securely in 
the range T^n = 6700-8100 K and are well-constrained, 
within the typical uncertainties of < 250 K, by the EW 
ratio measurement. One star (KIC 3230227) has EW ra¬ 
tios that indicate higher temperatures than the ranges 
covered by the MARCS models, so we have estimated 
the temperature by extrapolating (dashed line) the 7500- 
7750 K main-sequence points. Tables 5 (line 10) and 6 
(line 11 ) list our best-fit temperatures and uncertainties 
along with the temperatures estimated by TH12. Our 
temperature measurements typically he within 1 - 2(7 of 
those reported by TH12. KIC II494I30 is the most dis- 
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Figure 2. Effective temperature versus ratio of equivalent widths 
EWHa/EWsiii ■ Gray tracks show the measured equivalent width 
ratios for MARCS model atmospheres as a function of tempera¬ 
ture for five different surface gravities, as labeled from log 5- =3.0 
to log^ =5.0. The solid black line shows the nominal main se¬ 
quence, while the dashed line shows a quartic extrapolation of the 
7500-7750 K points used to estimate higher temperatures. Verti¬ 
cal lines mark the measured EW ratio of each star in our sample 
and provide a measure of temperature at various surface gravities. 
The temperature of each heartbeat star is estimated using a linear 
interpolation along the main-sequence track, denoted by the circle 
and error bars. (A color version of this figure is available in the 
online journal.) 

crepant; TH12 report a temperature of 6750 K, while 
we find a temperature of closer to 7500 K. Our EW ratio 
temperature estimate for KIC 3230227 is particularly un¬ 
certain because of the double-lined nature of the system 
(blended lines) and the fact that KIC 3230227’s tem¬ 
perature lies in the extrapolated regime. Consequently, 
we prefer a different method to estimate temperature as 
discussed in section 4.7.1. 

4. ORBITAL SOLUTIONS 

We used the power spectrum generated from a dis¬ 
crete Fourier transform of the radial velocity data to 
make an initial estimate of the orbital period for each 
system. Kiminki et al. (2012) provide an explanation 
of this method applied to massive O and B type stars. 
The power spectra for heartbeat stars are more complex 
than those for more “normal” binary systems that do 
not exhibit tidally induced pulsations. Power spectra for 
heartbeat stars show several strong peaks at periodici¬ 
ties of a few days to a few tens of days. The strongest 
peak is typically >50% higher than other peaks and, in 
all cases, correctly selects the orbital period of the sys¬ 
tem. These periods agree well with the more precise pe¬ 
riods inferred by TH12 from the the Kepler light curves, 
lending confidence to this independent determination of 
the orbital period. Other strong peaks (on the order 
of 2-4 days) likely correspond to the harmonic modes 
of the primary star induced by the close passage of the 
secondary star, as argued by TH12. We used the Bi¬ 
nary Star Combined Solution Package (Gudehus 2001) 
to solve for the orbital parameters and uncertainties, 
including a more refined period. Table 5 displays our 
best-fit results and the fitted parameters from TH12 for 
each SBl, listing the period (P) in days, the eccentricity 
of the system (e), the angle from node to periastron in 
degrees (w), the systemic velocity in km s“^ ( 7 ), the 


epoch of periastron (To), the projected velocity semi¬ 
amplitude for the primary in km s“^(A%), the projected 
semi-major axis of the system in AU (asinf), the mass 
function (M 2 sini)^/(Mi-|-M 2 )^ in Mq, the estimated 
inclination range in degrees (i), the estimated effective 
temperature of the primary (Peff), the estimated spec¬ 
tral type and adopted mass of the primary (Mi), the 
inferred secondary mass at the maximum estimated in¬ 
clination (M 2 ), the minimum mass ratio (g; calculated 
as M 2 at divided by Mi), and the rms of the residu¬ 
als of the best-fitting solution to the orbital elements in 
km s“^ (rmsi). The first eight rows of the table are di¬ 
rectly fitted with the Gudehus (2001) package while the 
next seven are estimated using methods described above 
and below. Table 6 lists the best-fit orbital elements 
and the TH12 ht for the SB2 systems, with the period 
(P) in days, the eccentricity of the system (e), the an¬ 
gle from node to periastron in degrees (w), the systemic 
velocity in km s“^ ( 7 ), the epoch of periastron (Tq), the 
projected velocity semi-amplitude for the primary (A%) 
and secondary {K 2 ) in km s“^, the projected semi-major 
axis of the system in AU (asini), the mass ratio q, the 
sum of the masses (Mi -|- M 2 ) sin^ i (M©), the estimated 
inclination range in degrees (i), the estimated effective 
temperature of the primary (Tett), the estimated spectral 
type and adopted mass of the primary (Mi) and sec¬ 
ondary (M 2 ), and the rms velocity of the residuals of the 
best-fitting solution to the primary and secondary (rmsi 
and rms 2 , respectively). 

4.1. Inclinations 

While SBl radial velocities alone cannot constrain sys¬ 
tem inclination, we can put limits on the inclinations of 
these systems and consequently place limits on the stellar 
masses. These masses and inclinations depend upon the 
adopted primary mass, an assumed stellar mass-radius 
relationship, and the accuracy of our measured orbital 
elements. There are two key pieces of information that 
allow inclination constraints: first, the mass of the sec¬ 
ondary cannot be larger than the mass of the primary, 
and second, these systems have been observed with Ke¬ 
pler to have either eclipsing or non-eclipsing behavior. 
The first criterion provides an upper limit on secondary 
mass and lower limit on the inclination. If a system is 
eclipsing, we can put a more stringent lower limit on the 
inclination by calculating the impact parameter (the pro¬ 
jected distance at conjunction relative to the observer), 
given by 

^ a cos f / 1 — \ 

Aprimary V 1 + ^ Siu W / 

(Winn 2010). Here, 6=1 indicates where one would 
begin to observe a grazing eclipse. If a system is not 
eclipsing, we estimate an upper limit for the inclination 
and a lower limit on the mass based on the inclination 
at which we would expect the system to show an eclipse. 
We calculate the impact parameter by adopting a main 
sequence dwarf mass consistent with our observed tem¬ 
perature and then applying a stellar mass-radius rela¬ 
tionship. We examined mass-radius relations of both 
Demircan & Kahraman (1991) and Maikov (2007), and 
there is excellent agreement in the calculated masses and 
inclinations between the two prescriptions (less than 1.5° 
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Table 5 

Orbital Elements of Single-Lined Binaries 


Element 

KIC 4248941 

KIC 8719324 

KIC 9899216 

KIC 11494130 


This work 

TH12 

This work 

TH12 

This work 

TH12 

This work 

TH12 

P (days) 

8.645 ± 0.002 

8.645 ± 0.001 

10.235 ± 0.005 

10.233 ± 0.001 

11.987 ±0.018 

10.916 ±0.001 

18.973 ± 0.008 

18.956 ± 0.001 

e 

0.34 ±0.04 

0.42 ±0.01 

0.64 ± 0.05 

0.60 ±0.01 

0.66 ± 0.24 

0.65 ±0.01 

0.49 ± 0.05 

0.62 ±0.01 

UJ (deg) 

308 ±5 

310 ±4 

330 ± 7 

343 ± 1 

35 ± 16 

81 ± 1 

70 ±4 

62 ± 1 

7 (km s“^) 

4.4 ±1.1 


-11.7± 1.9 


-1.0 ± 1.0 


-8.6 ±0.7 


To (HJD-2,450,000) 

6095.3 ±0.1 


6292.9 ±0.1 


6296.7 ±0.4 


6289.7 ±0.1 


Ki (km s“^) 

40.3 ±1.8 


43.5 ±3.1 


8.5 ±4.6 


27.0 ±1.0 


a sin 2 (AU) 

~0.03 


~0.03 


~0.01 


~0.04 


Mass function (M0) 

0.051 ±0.006 


0.065 ±0.013 


0.0003 ± 0.0002 


0.026 ± 0.002 


i (degrees) 

31-64 

68 ±6 

> 37 

76 ± 1 

6-15 

21 ± 1 

23-67 

48 ± 1 

Teff (K) 

6822 ± 107 

6750 

7597 ± 236 

7750 

7355 ± 201 

7500 

7456 ± 179 

6750 

Estimated S. C.i 

F3V 


A8V 


FOV 


A9V 


Ml (adopted; Mq) 

1.5 


1.7 


1.6 


1.7 


M 2 at iniax (Mq) 

0.7 


0.7 


0.5 


0.6 


Minimum q 

0.46 


0.43 


0.28 


0.33 


rmsi (km s“^) 

3.33 


6.44 


2.50 


3.84 



Table 6 

Orbital Elements of Double-Lined Binaries 


Element 

KIC 3230227 

This work TH12 

KIC 3749404 

This work TH12 

P (days) 

7.051 ±0.001 

7.047 ±0.001 

20.316 ±0.010 

20.307 ± 0.001 

e 

0.60 ± 0.04 

0.59 ± 0.01 

0.66 ± 0.03 

0.64 ± 0.01 

UJ (deg) 

293 ±4 

292 ± 1.2 

308 ±3 

302 ± 1 

7 (km s“^) 

-15.7 ± 1.7 


-13.9 ±0.8 


To (HJD-2,450,000) 

6311.76 ±0.03 


6656.85 ±0.10 


Ki (km s“^) 

98.5 ± 5.4 


50.0 ±2.1 


K 2 (km s“l) 

104.9 ± 6.1 


78.1 ±5.7 


a sin 2 (AU) 

~0.05 


~0.07 


g 

0.95 ±0.05 


0.64 ±0.09 


(Ml + M2) sin^ i (Mq) 

3.1 ±0.3 


1.9 ±0.3 


i (degrees) 

66-71 

43 ± 1 

51-65 

38 ± 1 

Teff (K) 

'^8000 

8750 

7231 ± 166 

7500 

Estimated S. C.i 

A6V 


FOV 


Ml (adopted; Mq) 

2.0 


1.6 


Estimated S. C.2 

A6V 


GIV 


M 2 (Mq) 

2.0 


1.0 


rmsi (km s~^) 

12.11 


2.01 


rms2 (km s~^) 

11.14 


0.85 



in inclination and 0.04 Mq in mass). Tables 5 and 6 and 
the sections below contain the individual results for each 
star. 

4.2. KIC 4248941 

We observed the SBl KIC 4248941 on 21 nights from 
2012 March through 2014 May. Our temperature esti¬ 
mate of T = 6822 ± 107 K agrees with TH12’s tempera¬ 
ture of 6750 K. We find a best-fit period of P = 8.645 ± 
0.002 days, which corresponds exactly to the TH12 pe¬ 
riod. The eccentricity is fit as e = 0.34 ± 0.04, which is 
similar to TH12’s reported eccentricity of 0.42. We also 
find a systemic velocity of 7 = 4.4± 1.1 km s“^ and a ve¬ 
locity amplitude of Ki = 40.3± 1.8 km s“^. This system 
has a projected semi-major axis of asini = 0.03 AU. 
Figure 3 displays the best-fit orbital solution, phased 
velocity curve, and phased light curve for this system. 
Black circles denote the radial velocity measurements ob¬ 
tained with WIRO, and the solid black line shows our 
best-fit solution. The highly structured curve plotted 
in the background shows the phased light curve taken 
from the publicly released simple aperture photometry 
(SAP) data in the Kepler MAST Archive, which was 


then phased and normalized with PyKE (Still & Barclay 
2012). The overplotted dashed curve shows the Kumar et 
al. (1995) light curve model generated using our orbital 
parameters and the TH12 inclination.^ The horizontal 
axis shows the orbital phase; the left vertical axis marks 
the velocity of the measured components in km s“^, while 
the right vertical axis shows the normalized flux of the 
light curve. Note that the point of maximum velocity 

^ The model of photometric tidal variations was taken from equa¬ 
tion 44 in Kumar et al. (1995), presented here as written in equa¬ 
tion 1 of TH12: 

1 - isirP - uj) ^ 

f ~ ' im/af 

Here, ^ is the change in flux of the light curve, 5 is a scaling 
factor to match the amplitude of the normalized light curve, i is 
the system inclination (for our purposes, taken from TH12), (f){t) 
is the true anomaly (calculated by first computing the eccentric 
anomaly, which itself is dependent upon the eccentricity), uj is the 
angle from node to periastron, R{t) is the time-variable separation 
between the two stars (dependent upon a, e, and the eccentric 
anomaly), a is the semi-major axis, and C is an offset to represent 
the zero-point of the normalized light curve. 
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Figure 3. The radial velocity data and best-fit solution for the 
8.645 day period SBl system KIC 4248941. Circles denote the 
primary star’s velocity. The left vertical axis shows Heliocentric 
radial velocity, and the horizontal axis shows the orbital phase. 
The highly structured curve plotted in the background is the Kepler 
light curve with the normalized flux on the right vertical axis. A 
model light curve created with the orbital elements from the data in 
this paper and the TH12 inclination is shown in the dashed curve. 
(A color version of this figure is available in the online journal.) 



Figure 4. The radial velocity data and best-fit solution for the 
10.235 day period SBl system KIC 8719324. This plot uses the 
same conventions as Figure 3. (A color version of this figure is 
available in the online journal.) 

amplitude (at phase=l) coincides with largest variation 
in the light curve. This is consistent with the scenario 
proposed by TH12 in which the primary star undergoes 
large tidal distortions near the time of periastron. 

To estimate the mass of the secondary, we assume a 
primary spectral type of F3V (1.5 Mq) and no eclipses. 
We find that the inclination must lie in the range 31- 
64°, corresponding to a secondary masse of 1.5-0.7 Mq. 
Thus, the minimum mass ratio is g = 0.46. TH12’s in¬ 
clination of 68±5.5° is consistent with our range. 

4.3. KIC 8719324 

KIC 8719324, an SBl system, was observed on 24 
nights from 2012 March through 2014 May. We find a 
temperature of T = 7597 ± 236 K, which is similar to 
TH12’s value of 7750 K. This system shows a grazing 
eclipse signature in the light curve as reported in TH12, 
although the light curve is dominated by the signature of 
the KIC 8719324 tidal modes. We find a period of P = 


Figure 5. The radial velocity data and best-fit solution for the 
11.987 day period SBl system KIC 9899216. This plot uses the 
same conventions as Figure 3. (A color version of this figure is 
available in the online journal.) 

10.235 ± 0.005 days, an eccentricity of e = 0.64 ± 0.05, a 
systemic velocity of 7 = —11.7 ± 1.9 km s“^, a velocity 
amplitude of Ki = 43.5 ± 3.1 km s“^, and a projected 
semi-major axis of asini = 0.03 AU. Our solution agrees 
well with the TH12 period of 10.233 days and eccentric¬ 
ity of 0.60. Figure 4 depicts the best-fit solution, data, 
and the Kepler light curve for KIC 8719324. 

We estimate the mass of the secondary to be a mini¬ 
mum of 0.7 Mq by assuming an A8V spectral type with 
mass 1.7 Mg for the primary star. By using the pres¬ 
ence of a grazing eclipse, we constrain the inclination to 
be greater than 37°, which limits the secondary mass to 
less than 1.5 Mq. This inclination range supports the 
TH12 value of 76° and provides a minimum mass ratio 
of g = 0.43 at i = 90°. 

4.4. KIC 9899216 

KIC 9899216, another SBl, has an effective temper¬ 
ature of T = 7355 ± 201 K, in agreement with TH12’s 
temperature of 7500 K. We have sixteen nights of ob¬ 
servation for this object from 2013 May to 2014 June. 
The spectroscopic period of this system is P = 11.987 ± 
0.018 days, which is somewhat larger than TH12’s pe¬ 
riod of 10.916 days. This discrepancy in periods likely 
arises because of the low velocity amplitude of Ki = 
8.5±4.6 km s“^, which is roughly twice the typical veloc¬ 
ity uncertainty on each measurement. We did attempt to 
fix our orbital solution at the period of TH12, but this re¬ 
sulted in a much poorer solution. Additional data will be 
required to reconcile the difference in periods. However, 
we do find the eccentricity of e = 0.66 ± 0.24 to be com¬ 
parable to TH12’s eccentricity of 0.65. We also measure 
the systemic velocity to be 7 = —1.0 ± 1.0 km s“^ and 
the projected semi-major axis to be a sin* = 0.01 AU. 
Figure 5 depicts the best-fit solution, data, and Kepler 
light curve for KIC 9899216. 

TH12 estimate that KIC 9899216 has an inclination 
of 21°. This is the lowest estimated inclination of our 
sample of heartbeat stars. Our inclination range of 6- 
15° corroborates the TH12 result. The minimum sec¬ 
ondary mass of 0.5 Mq, coupled with the estimation that 
KIC 9899216 is an FOV star with mass 1.6 Mq, gives a 
minimum mass ratio of g = 0.28. 
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Figure 6. The radial velocity data and best-fit solution for the 
18.973 day period SBl system KIC 11494130. This plot uses the 
same conventions as Figure 3. The discrepant data point at phase 
0.16 is not included in the fit. 

4.5. KIC 11494130 

We gathered 19 nights of data on KIC 11494130 from 
2012 June through 2014 May. We find a temperature of 
T = 7456 ± 179 K, which contrasts greatly with TH12’s 
temperature of 6750 K. TH12 fit this star with a logg of 
4.5; if we adopt the same value, our temperature would 
drop to near 7200 K. A log 5 of 5.0 yields a temperature 
much closer to the TH12 value, but is significantly larger 
than the canonical values for a main sequence star of 
this temperature. The best-fit orbital solution yields a 
period of P = 18.973±0.008 days, which is similar to the 
TH12 period of 18.956 days. The best-fit eccentricity is 
e = 0.49±0.05, somewhat less than the TH12 eccentricity 
of 0.62. However, our best-fit eccentricity is plausible 
if TH12 underestimated the width of the leading wing 
of the light curve pulsation, which is possible with the 
significant variability near periastron. We find a systemic 
velocity of 7 = — 8.6 ± 0.7 km s“^, a primary velocity 
amplitude of Ki = 27.0 ± 1.0 km s“^, and a projected 
semi-major axis of asini = 0.04 AU. Figure 6 displays 
the velocity data and best-fit solution, in addition to the 
Kepler light curve. The discrepant data point at phase 
0.16 is not included in the fit. 

Depending on the actual temperature of the star, the 
mass of the primary will differ. If we assume an A9V 
spectral type (7500 K and mass 1.7 Mq, from this work) 
for the primary, we find an inclination range of 23-67°. 
This leads to a minimum secondary mass of 0.6 Mq and 
a minimum mass ratio oi q = 0.33. Assuming instead 
an F5V spectral type for the primary (6700 K and mass 
1.4 Mq, from TH12), we find an inclination range of 
25-70°, giving a minimum secondary of mass 0.5 Mq. 
However, this leads to the same minimum mass ratio. 

4.6. KIC 3749404 

KIC 3749404 is an SB2 for which we have eighteen ob¬ 
servations between 2013 May to 2014 September. We 
find an effective temperature of T = 7231 ± 166 K, which 
agrees well with TH12’s reported temperature of 7500 K. 
The best-fit period of P = 20.316 ± 0.010 days matches 
TH12’s photometric period of 20.307 days. The eccen¬ 
tricity was best fit with a value of e = 0.66 ± 0.03, con¬ 
sistent with the TH12 estimate of 0.64. Our fit also pro¬ 
duced a systemic velocity of 7 = —13.9 ± 0.8 km s“^ 



Figure 7. Radial velocity data and best-fit solution for the 20.316 
day period SB2 system KIC 3749404. Circles denote the primary, 
while triangles show data for the secondary component. This figure 
follows the same conventions as Figure 3. (A color version of this 
figure is available in the online journal.) 

and a velocity amplitude of Ki = 50.0 ± 2.1 km s“^ for 
the primary and K 2 = 78.1 ± 5.7 km s“^ for the sec¬ 
ondary. This system has a projected semi-major axis 
of asinz = 0.07 AU. Hambleton et al. (2013) note that 
the light curve shows an apsidal advance over the four 
years of Kepler observations, which may lead to an un¬ 
derstanding of the binary evolutionary state. 

Line splitting in this system occurs only near perias¬ 
tron, so just two measurements of the secondary veloc¬ 
ity can be obtained with our spectra. These data are 
enough, though, to place good constraints on the mass 
of the secondary star. The mass ratio is measured to 
be <7 = 0.64 ± 0.09 and (Mi -|- M 2 )sin^z, the sum of 
the masses times the sine of the inclination cubed, is 
1.9 ± 0.3. If we adopt an FOV spectral type correspond¬ 
ing to a mass of 1.6 Mq for the primary, we find that the 
secondary star has a mass of 1.0 Mq. This corresponds 
roughly to an GIV spectral type. The inclination of the 
system can then be calculated to lie in the range 51-65°, 
significantly higher than TH12’s reported inclination of 
38°. 

If the primary is an FOV star with a GIV secondary, 
we would expect roughly 80% of the flux to come from 
the primary and 20% from the secondary. The redshifted 
component we see at periastron has an equivalent width 
in Si II A6371 that is about 20-25% of the total equivalent 
width we see in a blended spectrum. Thus, the secondary 
mass we derive above is consistent with the flux ratios 
seen in our spectra. 

4.7. KIC 3230227 

Our spectra reveal that KIG 3230227 is a double-lined 
spectroscopic binary. TH12 reported grazing eclipses in 
the light curve but did not note the double-lined behav¬ 
ior. Our observations include 23 nights from 2012 March 
through 2013 June. We find that the two components are 
of roughly equal temperature, estimated to be 8000 K 
(discussed below), while TH12 estimated a temperature 
for the system closer to 8750 K. 

For KIG 3230227 our best-fit solution provides a pe¬ 
riod of P = 7.051 ± 0.001 days. This agrees well with 
the TH12 value of 7.047 days. We find an eccentric¬ 
ity of e = 0.60 ± 0.04, in good agreement with TH12’s 
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Figure 8. The radial velocity data and best-fit solution for 
KIC 3230227, a 7.051 day period SB2 system. Circles denote the 
primary, while triangles are the data for the secondary component. 
This figure follows the same conventions as Figure 3. The grazing 
eclipse at phase=1.0 (which has a 6% depth in the full light curve) 
has been truncated 1% to allow the stellar pulsations to be seen. 
(A color version of this figure is available in the online journal.) 


value of 0.59. The velocity amplitude of the primary 
is Ki = 98.5 ± 5.4 km s“^, while the secondary has 
an amplitude of K2 = 104.9 ± 6.1 km s“^. The sim¬ 
ilarity of the velocity amplitudes indicates a mass ra¬ 
tio very close to unity. We find a systemic velocity of 
7 = —15.7±1.8 km s“^ and a projected semi-major axis 
of asinz = 0.05 AU. Figure 8 shows the best-fit solution 
and measured velocities. Circles denote data from the 
primary while triangles denote data from the secondary. 
There is a 6 % decrease in normalized flux at the pri¬ 
mary eclipse; the secondary eclipse is not seen owing to 
the highly eccentric orbit. The primary eclipse has been 
truncated at 1% of its observed depth in Figure 8 to al¬ 
low the stellar pulsations to be more easily seen. Our 
analytic light curve model (dashed line) does not include 
an eclipse and may have a small scaling error due to un¬ 
certainties in the depth of the tidal variation. 

Figure 9 shows a time series of the Si II A6347 and Ha 
lines ordered by orbital phase with the horizontal axis in 
km s“^ and the dotted line at 0 km s“^ showing the rest 
wavelength. Both lines display an obvious splitting near 
phase=0. This splitting is mimicked in other lines, such 
as the Si II A6371 feature. Equivalent widths of the sili¬ 
con and Ha features, as well as other lines not pictured, 
appear slightly smaller in the secondary component, but 
the spectral types must be similar because the spectral 
features are nearly identical in unblended spectra. 


4.7.1. Temperature 

Because of the difficulty in calculating temperatures 
for the individual components of KIC 3230227, we de¬ 
cided to combine model atmospheres in an attempt to 
reproduce the observed spectrum at periastron when the 
two stellar components are minimally blended. We used 
the PHOENIX model atmospheres of Husser et al. (2013) 
because the spectral resolution is very similar to that of 
our data and the temperature grid spans the range of 
A-type stars in increments of 200 K. We corrected the 
model spectra to air wavelengths and continuum nor¬ 
malized the flux in a manner similar to observed spectra. 
We shifted the blue component to —146.5 km s“^ and 



Figure 9. Profiles of the Si II A6347 feature and Ha in 
KIC 3230227 when ordered by phase (labeled under each spec¬ 
trum). The horizontal axis shows velocity centered at the rest 
wavelength of the feature, as denoted by the vertical dotted line. 
Both features show obvious splitting at phases near periastron pas¬ 
sage (0 and 1). 


the red component to 122.6 km s“^ to reflect the maxi¬ 
mum velocities in our observed spectra. We also added in 
an empirical Na I D profile modeled after the interstellar 
sodium absorption of a spectrum where the lines are fully 
blended. Finally, we multiplied our resultant spectrum 
by the atmospheric transmission spectrum from Hinkle et 
al. (2003). We find the atmospheric transmission neces¬ 
sary to include because the H 2 O feature near the sodium 
doublet and the O 2 feature close to the silicon lines, to 
name two, affect the continuum and depth of the ob¬ 
served features. Although the atmospheric transmission 
at Kitt Peak measured by Hinkle et al. (2003) may not 
exactly match that at WIRO, this zeroth order correc¬ 
tion is sufficient to show the probable magnitude of the 
atmospheric features in the spectra. 

Figure 10 shows the best-fit synthetic model spectrum 
(bold and dotted line; colored red in the electronic edi¬ 
tion) overplotted on the observed spectrum (thin black 
line). The lower panel shows the full spectral range of 
our data; the upper three panels show magnihed views 
of important features: Na I on the left, Si II in the mid¬ 
dle, and Ha on the right. The rest wavelength of each 
line has been marked as “0 km s“^”, and the velocities of 
±150 km s“^ are also labeled. In the Na I panel, the red- 
shifted component of the A5889 line and the blueshifted 
component of the A5895 line overlap. From our best- 
fit synthetic spectrum, we estimate that the stars have 
very similar spectral types near A6-7V. Both model at¬ 
mospheres shown have a temperature of 8000 K. The 
redshifted secondary component has a surface gravity of 
logg = 3.5 and rotational broadening of 75 km s“^, while 
the blueshifted primary component has a surface gravity 
of logg = 4.0 and rotational broadening of 30 km s“^. 
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Overall, the figure shows good agreement between the 
synthetic spectrum and the observed spectrum, espe¬ 
cially in regions like Ha and Si II. However, although 
the ratios of the sodium lines in the model atmosphere 
match those in the observed spectrum, we have not been 
able to replicate the absolute depth despite having varied 
parameters such as rotational broadening, spectral type, 
and surface gravity. Increasing the contribution of the 
interstellar Na does not improve the fit because the in¬ 
terstellar lines fall at zero velocity and are much weaker 
than the stellar lines. Adopting the hotter effective tem¬ 
perature found by TH12 only serves to exacerbate the 
discrepancy between the synthetic and observed spectra. 
Increasing the metallicity to [-1-0.5] and [-1-1.0] resulted in 
a slightly better fit to the sodium lines but a much poorer 
fit to the silicon and other metal lines, in the sense that 
the metal-rich models greatly overpredicted the strength 
of metal lines throughout the spectrum. Thus, it is pos¬ 
sible that these stars are more metal rich than solar or 
have peculiar metal abundances, but the difference is not 
reproducible with the available PHOENIX spectra. 

4.7.2. Inclination 

Knowing the orbital parameters for the stars of 
KIC 3230227, we calculate that the sum of the masses 
times sin^ i is 3.1 ± 0.3 Mq and that the mass ratio is 
q = M2/M1 = 0.94 ± 0.05. Adopting TH12’s inclina¬ 
tion of 43°, the sum of the masses would be ^10 Mq. 
This would correspond to two late-type B stars, in strong 
disagreement with the fitted temperatures and spectral 
types. An inclination of 66-71° is required to yield a sum 
of ^4 Mq appropriate to two mid-type A stars. Inclina¬ 
tions as low as 55° and as high as 80° may be allowed 
owing to the uncertainties on the components’ spectral 
types. TH12 note that their calculated inclination for 
KIC 3230227 seems low for an eclipsing system and that 
it may be underestimated. The most likely reason, ac¬ 
cording to their analysis, is because the model from Ku¬ 
mar et al. (1995) does not include stellar irradiation. 

5. DISCUSSION AND SUMMARY 

We have presented spectroscopic orbital solutions for 
six heartbeat stars, further supporting the hypotheses of 
Welsh et al. (2011), Thompson et al. (2012), and others 
that this class of stars with highly variable light curves 
are indeed high-eccentricity binary systems. All six 
systems have best-fit eccentricities > 0.34 (most above 
0.5) and periods of 7-20 days. The good agreement 
between our spectroscopic orbital solutions and the so¬ 
lutions derived from Kepler light curves in TH12 attests 
to the applicability of the Kumar et al. (1995) model of 
static modes of gravitationally-induced tidal variations 
to heartbeat stars orbital element extraction. 

Each system in our sample contains an A- or F-type 
primary component. The six systems, of which two are 
SB2s and four are SBls, have minimum mass ratios 
q ~ 0.3-1.0, making the secondaries probable M-dwarf 
stars or earlier. KIC 3230227 is a double-lined system 
containing two mid-A stars with a mass ratio very close 
to unity. This system is reminiscent of KOT54, the pro¬ 
totypical heartbeat system (Welsh et al. 2011), which 
consists of two massive A-type stars in a 41-day orbit. 
Our other SB2, KIC 3749404, exhibits similarities to 
KIC 3230227, such as an eccentricity on the higher end 


of our sample, but has a smaller primary mass, longer 
period, and smaller mass ratio. The two SB2s have the 
highest mass ratios of our sample, and both SB2s have in¬ 
clinations higher than reported in TH12. This may be an 
indication of stellar irradiation or other effects creating 
signatures in the light curve that are not included in the 
Kumar et al. (1995) model. These limited data suggest 
that photometrically derived orbital elements systemat¬ 
ically underestimate the inclination for high mass ratio 
systems, although more data would be required to con¬ 
firm this observation. The heartbeat stars herein exhibit 
mass ratios q > 0.3, suggesting that relatively massive 
companions may be required to excite the tidal and pul- 
sational modes that are the hallmark of these intriguing 
systems. However, a wide range of companions are al¬ 
lowed. On the basis of this limited sample size, there 
does not appear to be evidence for a correlation between 
primary mass and mass ratio or between orbital elements 
and mass ratio. These spectroscopic observations pro¬ 
vide fundamental data on the masses and mass ratios of 
heartbeat stars that can be used to further understand 
their complex photometric variability, orbital dynamics, 
and binary evolution. 
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Figure 10. The best-fit spectrum for KIC 3230227. The bold line depicts the composite model spectrum, consisting of a blueshifted 
model stellar atmosphere with a temperature of 8000 K and a surface gravity of \ogg = 3.5, a redshifted model atmosphere with T=8000 K 
and logg = 4.0, an empirical model of the interstellar sodium doublet Na I D AA5889,5995, and the atmospheric transmission spectrum of 
Hinkle et al. (2003) smoothed to the resolution of the data. The thin black line shows a spectrum taken with WIRO when the components 
are at maximum separation. The top three panes show magnified versions of the Na I D doublet, the Si II lines, and Hu, respectively. 
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feature overlap. The bottom panel shows the full spectrum with the above three features marked. We also label two atmospheric regions 
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available in the online journal.) 
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